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Abstract 

We derive a compact expression for the three-point MHV form factors of half-BPS op- 
erators in M = 4 super Yang-Mills at two loops. The main tools of our calculation are 
generalised unitarity applied at the form factor level, and the compact expressions for 
supersymmetric tree-level form factors and amplitudes entering the cuts. We confirm 
that infrared divergences exponentiate as expected, and that collinear factorisation is en- 
tirely captured by an ABDK/BDS ansatz. Next, we construct the two- loop remainder 
function obtained by subtracting this ansatz from the full two-loop form factor and com- 
pute it numerically. Using symbology, combined with various physical constraints and 
symmetries, we find a unique solution for its symbol. With this input we construct a 
remarkably compact analytic expression for the remainder function, which contains only 
classical polylogarithms, and compare it to our numerical results. Furthermore, we make 
the surprising observation that our remainder is equal to the maximally transcendental 
piece of a closely related two- loop amplitude in QCD. 
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1 Introduction 



The last years have witnessed dramatic progress in our understanding of seemingly unre- 
lated physical quantities in M = 4 supersymmetric Yang-Mills theory (SYM), in particular 
scattering amplitudes, Wilson loops and correlations functions. One important lesson is 
that these objects are actually related in unexpected ways, which has led to major con- 
ceptual insights, vastly improved techniques to calculate them and the discovery of new 
structures and symmetries. 
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In this context it is natural to consider simple, interesting generalisations of the quanti- 
ties mentioned above, such as form factors. In a gauge theory, one considers the overlap of 
a state created by a gauge- invariant operator, 0(x) |0), with a multiparticle state (1 • • • n|. 
The form factor is then defined as 

/n 
d*xe~ iqx (l---n\O(x)\0) = 5 {4) (q - F(l, . . . , n) , (1.1) 

i=i 

where F(l,...,n) := (1 • • • n|O(0)|0). 

For n = 2, the latter is the celebrated Sudakov form factor, which is a famous example 
of exponentiation - a phenomenon discovered more recently in scattering amplitudes of 
M = 4 SYM [JJ[2], and is intimately linked to scattering amplitudes, whose universal 
infrared divergences are captured by Sudakov form factors. An n-point form factor can 
also be viewed as describing the decay of a state created by the gauge-invariant operator 
O into an n-particle state. For n = 2 the kinematic dependence is trivial, and the form 
factor can be expressed in terms of the cusp and collinear anomalous dimensions. In 
planar M = 4 SYM, the former is known to all orders in the gauge coupling [3]. 

In this paper we will mainly be concerned with form factors with n > 2, which have 
highly non-trivial kinematic dependence. While form factors share several properties with 
scattering amplitudes, such as soft and collinear factorisation, they also display interesting 
and important differences. This is due to the insertion of a gauge-invariant operator, which 
carries momentum that can be injected at any location in Feynman diagrams, and is thus 
not restricted by the colour ordering of the external, on-shell particles. In practice this 
leads to a more complicated colour management and to the appearance of non-planar 
loop integral topologies, even if we restrict ourselves to the planar limit, as we do in this 
paper. Furthermore, this also implies that the dual conformal (super) symmetry of planar 
amplitudes in M = 4 SYM is not a property of form factors. 

Form factors at strong coupling in M = 4 SYM were recently studied in [HE]. The 
particular choice of the operator considered in those papers is immaterial, as long as 
its conformal dimension is small compared to VX Specifically, in [5] form factors were 
analysed in great detail and the explicit expression of the four-point form factor was 
determined in the case of (1 + l)-dimensional kinematics. This is the first example of a 
nontrivial form factor remainder at strong coupling. In [6l[7] we have begun a systematic 
investigation of form factors of bilinear half-BPS operators at weak coupling, following 
a pioneering paper of Van Neerven [8jJ3 In particular, in [6] we determined the explicit 
one-loop expressions for these form factors for an arbitrary number of external particles 
with an MHV helicity configuration. 

In this paper we will focus on the two-loop calculation of the three-point form factor. 
In QCD, similar quantities have been calculated at one [H] and, more recently, two loops 
P2J[13] using Feynman diagrams. These form factors are phenomenologically important 
since they are related to the scattering of e + e~ — > 3 jets and H — > 3 jets. Here we will use 
instead (generalised) unitarity [141417] applied directly at the level of the form factor [6] 
in order to determine the two-loop, three-point form factor from its cuts. 

1 A parallel investigation was carried out in [9lll0j. 
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In Section 2 we will set the stage by presenting the details of a simpler calculation, 
namely that of the Sudakov form factor at two loops. Section 3 describes the calculation 
of our three-point form factor at two loops from generalised cuts. We obtain a compact 
result, presented in (13. lip , expressed in terms of planar and non-planar two-loop integral 
functions. Most (but not all) of the analytic expressions for these functions have been 
derived in [T8H20] . For this reason, we have resorted to numerical methods in order to 
collect data about this form factor which can then be compared to analytic expressions. 
The main tool here is the Mellin-Barnes representation of higher-loop integrals (see for 
example [2T]), and their evaluation using the powerful numerical algorithm of |22j. 

In order to present the two-loop form factor in an efficient way, we construct in Section 
4 a remainder function very much in the way as for the case of scattering amplitudes. At 
two loops, this remainder is 

^ ■= Si 2 ^) " \(G { n 1) (e)) 2 -f^(e)g^(2e) - C7< 2 > + 0(e) , (1.2) 

where is the colour-stripped L-loop form factor divided by its tree-level expression, 
is an n-independent constant and /^(e) contains the cusp and collinear anomalous 
dimensions. One important reason to introduce this remainder is that, as we show ex- 
plicitly in Section 4, it has the correct collinear limits. More precisely, if is chosen 
appropriately then TZ^ = and in a simple collinear limit TZn — > TZn-i- m other words, 
the ABDK/BDK ansatz for the form factor captures correctly the collinear behaviour of 
the two-loop form factor. The remaining task is then to determine this remainder function 
which is a function of scale invariant ratios of Mandelstam variables only. 

Our main weapon to attack this remainder analytically is the use of symbols |23j . 
This is a very powerful concept that allows to re-express complicated identities between 
polylogarithmic functions in terms of simple algebraic identities. It was used in [23] to 
rewrite the Wilson loop six-point remainder function found in |24j in a very compact 
form. Because of the amplitude/ Wilson loop duality [25H27] . this remainder is also equal 
to the six-point MHV amplitude remainder at two loops, as was confirmed numerically 
in [28l[29]. Recent interesting applications of the symbol to the construction of various 
remainder functions in M = 4 SYM were presented in [301435] . 

As we will discuss in Section 4, we can impose several physical constraints, such as 
collinear limits and various requirements on the first, second and last entry of the symbol, 
directly at the level of the symbol. In the present instance, we find that the various 
physical properties are so stringent that the symbol of the two-loop remainder function 
at three points is determined uniquely; we present its explicit expression in (14.281) . An 
important point we would like to anticipate is that the symbol of the three-point form 
factor at two loops depends on a very restricted set of variables, as we discuss in detail in 
Section 4. Furthermore, it obeys a particular symmetry condition [23] which guarantees 
that the corresponding function can be written in terms of classical polylogarithms up to 
degree four only. Using this information, we will construct a surprisingly simple analytic 
expression for the complete remainder function, shown in (I4.32p . We note that beyond- 
the-symbol ambiguities can easily be fixed by imposing correct behaviour in collinear 
limits. Within errors, our result is in good agreement with our numerical results obtained 
for various kinematic points. 



3 



We also observe an interesting relation between the symbol of the three-point form 
factor remainder and a piece of the six-point MHV amplitude remainder of [23] upon 
appropriate identification of the kinematic variables. This is reminiscent of a result found 
in [5], where the four-point form factor in (1 + l)-dimensional kinematics was evaluated 
at strong coupling and expressed in terms of the eight-point MHV amplitude remainder 
in (1 + l)-dimensional kinematics. 

We will conclude by establishing an intriguing relation between our remainder and the 
finite remainder function (defined with a prescription introduced by Catani [36]) of the 
two-loop amplitude H — > ggg in QCD, a quantity that was recently presented in [T2|[T3] 
in the large top mass limit. Interestingly, if we translate the latter into a ABDK/BDS 
remainder function and extract the terms of highest transcendentality, which is four at 
two loops, we find a function that has the same symbol as the remainder of the three- 
point form factor calculated in this paper. In fact, this observation provides us with an 
alternative way to find the analytic M = 4 SYM remainder, although our expression is 
much more compact - specifically, the analytic formula obtained pursuing this alternative 
strategy contains also two-dimensional harmonic polylogarithms (2dHPL's) [TUPD] which, 
as we mentioned, are not present in our expression (I4.32p . Remaining beyond-the-symbol 
ambiguities can again be fixed using collinear limits and symmetry. We have checked 
by explicit comparison that these two representations of the remainder are in complete 
agreement. 

To our knowledge, this is the first occurrence of the principle of maximal transcenden- 
tality [37], which relates the maximally transcendental part of a quantity in QCD to the 
same quantity in planar M = 4 SYM, in an observable with non-trivial kinematic depen- 
dence. It would be interesting to find more examples of this, and understand whether this 
observation can be used to simplify the calculation and expressions of other scattering 
amplitudes in QCD. 



2 Warm up: the Sudakov form factor at two loops 

In this section we present a very simple, unitarity-based derivation of the two-loop Sudakov 
form factor F(q 2 ) := (0 12 (pi)0i 2 (p2)|Tr(0i 2 0i2) (0)|0), where q := pi + p 2 - This quantity 
was calculated at one and two loops in = 4 SYM in [S] using unitarity applied to 
Feynman diagrams [42j, see also [6_ll38l[39], and recently at three loops in [3D]. Here 
we will use an approach based on unitarity directly applied at the level of the form 
factor. This approach was introduced in [6J and is a straightforward generalisation of the 
unitarity-based method of [I3|[To]. 

There is only one kinematic channel, the q 2 channel, hence it is sufficient to consider 
two-particle cuts, and lift the cut integral to a full .D- dimensional integral, as in |45j . 
There are two distinct configurations contributing to such a cut: in the first one, a tree- 
level, two-point form factor and a one-loop amplitude enter the cut, whereas in the second 
one we have a one-loop two-point form factor, and a tree-level four-point amplitude. Let 
us focus on the first possibility. Furthermore, for the sake of illustration, we will perform 
the calculation in two different ways: firstly, using the conventional representation of the 



4 





(a) (b) 

Figure 1: The two two-particle cuts contributing to the Sudakov form factor. 

one-loop amplitude entering the cut in terms of fundamental colour generators, and then 
using its representation in terms of adjoint generators [4"6] . 



2.1 Calculation with fundamental generators 

The first quantity entering the cut is the four-point one-loop amplitude, whose expression 
in terms of fundamental colour generators is given in (1A.ll) with n = 4. This expression 
has to be convoluted with a tree- level two-point form factor, simply given by <5 ai i ai 2, see 
Figure [T]a. We focus first on the contribution from the planar amplitude, i.e. the first 
line in flA.lj) . There are six possible permutations to consider, which give rise to the 
single-trace structured 



Tr(l, 2, h, l 2 ), Tr(l, 2, l 2 , h), Tr(l, l u l 2 , 2), Tr(l, l 2 , h, 2), 
Tr(l,h,2,l 2 ), Tr(l,l 2 ,2,h). 



(2.1; 



When contracting with the tree form factor {(f>^ 2 (/i)0 12 2 (p2)|Tr(0i 2 0i 2 ) (0) 1 0) = 5 ai i a! 2, 
and using 5 ab T^ m = 5 im 5ji and b ah T^T h jm = NS im we see that only the first line of f l2~Tj) 
is leading in colour. The four traces give an identical result, namely ATTr(aia 2 ) = N 2 5 aiCl2 . 
The contribution to the cut from the planar one-loop amplitude is then 



N 6 



2 jra 1 a2 



43(1, 2, M 2 ) + Ami,2,l 2 ,h) + AW(2,l,h,h) + At.{(2,l,l 2 ,h 



iW 



(!)/ 



(2.2) 



We now consider the contribution from the non-planar part of the one-loop amplitude 
(see section A for details). The COP for {a} = {2, 1} and {/?} = {3,4} (corresponding 
to c = 3 in (ICTl ) are 

(2,1,3,4), (2,3,1,4), (3,2,1,4), (1,2,3,4), (1,3,2,4), (3,1,2,4), (2.3) 

and the non-planar one-loop piece (second line of (lA.ip ) is 

A$l(l, ... 4) = Tr(l, 2)Tr(k, l 2 )A^(l, 2, h, l 2 ) + Tr(l, Z 1 )Tr(2, 1 2 )A®(1, k, 2, l 2 ) 

+ Tr(l,/ 2 )Tr(2,/ 1 )Ag(l,Z 2 ,2,/ 1 ) , (2.4) 

2 To keep the notation simple we define Tr(l, 2, l lr --) := Tr(r ai T a2 T 0! i • • • ) etc. 
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with 

A<g(l,2,M 2 ) = Ai;j(2,l,i 1 ,i 2 )+Ag(2,i 1 ,l,Z 2 )+4!J(i 1 ,2,l,Z 2 ) 

+ Ag(l, 2, / 1; / 2 ) + 4*1(1, / 1; 2, Z 2 ) + Ag(/ 1; 1, 2, Z 2 ) . (2.5) 

Contracting with the tree-level form factor, we see that the leading structures in colour 
are of the form <5 a 'i a, 2Tr(/i/ 2 )Tr(12) = N 2 5 aia2 . Collecting terms, we obtain a two-loop 
cut integrand 

F {2)aia2 (l%^ cut = 2N 2 5 a ^ JdLlPS(h,l 2 ;q) [4A$(1, 2, l u l 2 ) + A%\(1, h, 2, l 2 ) , 

(2.6) 

where the Lorentz invariant phase space measure in dimensional regularisation is 

dLIPS(Zi, h] q) := d D h d D l 2 8 + {ll)8 + {ll)8 D {h + l 2 + q) , (2.7) 

with D = 4 — 2e. The one-loop component amplitudes appearing in (12. 6p are given, to all 
orders in the dimensional regularisation parameter e by 

Ag(l, 2, h, l 2 ) = A® (01 2 (P1), 01 2 (P 2 ), 034(/l), Mh)) ^(l, 2, / l5 Z 2 ) 

A$(l,h,2,l 2 ) = A( o )(0 12 ( Pl ),034(/i),0i2(p 2 ),034(/ 2 ))F Om (l,/i,2,/ 2 ) (2.8) 

where F 0m (l, . . . , 4) = si 2 s 2 3/ 0m (l, . . . , 4) is a zero-mass box function. The tree ampli- 
tudes A( o )(0i 2 (pi),0i 2 (p 2 ),034(/i),034(/ 2 )) and A(°)(0i 2 (p 1 ),034(/ 1 ),0i 2 (p 2 ),03 4 (/ 2 )) can 
be extracted from Nair's superamplitude 



Amy := 9 n - 2 (2vr) 4 ^ ( £ A, A,) ( £ \Vi) f[ y-^y , (2.9) 

i=l ' ' 



where A n+ i = Ai, givin 



i=l i=l 



^ (O) (012(Pl),012(p 2 ),034(/l),034(/ 2 )) = , (2-10) 

^ (O) (^bl)^34(/l),01 2 (p 2 ),034(/ 2 )) = n^)^ 1 , 2 ]], n ■ (2.H) 



(1Zi)(Zi2)(2Z 2 )(Z 2 1) 



On the cut one easily finds 

(12) (hh) q 2 



(2-12) 



(2/ 1 )(/ 2 i) (i 2 + Pl y ' 

Using this relation we easily arrive at 

F^ ab (q 2 ) = N 2 5 ah F^ 2 \q 2 ), (2.13) 

with 

F^ 2 \q 2 ) = 2[4LT(g 2 ,e) + CT(g 2 ,e)l , (2.14) 



3 Note that s 12 s 2h A^ (0i 2 (pi), </>i2(p2),<t>3i(h), 034 (fe)) = s Ul si l2 A^ (0i 2 (pi), <t>3i{h), ^lafe), 034(k)), 
which simplihes the cut algebra. 
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where the two-loop ladder and crossed triangle, LT(q 2 , e), CT(g 2 , e) respectively, are given 

by mHniafl 



where 



LT(g 2 ,e) 



CT(g 2 ,e) 



G(x,y) 
G 3 (x,y,z) 



2\-2e 2-Ye 

-q ) e ' 



-G(2,2)G 3 (2 + e,l,l) 



-G(2,l)[-G 3 (2,l,l + e) + G 3 (l,l,l) 



(-Q 2 ) 



2\-2e 



1 5tt 2 29 , 3 4 ^, N 



,2\-2e 



1 7T 2 83 59 

e 1 ~ e 2 " ~ 3e 3 ~ 120 



vr 4 + 0(e) 



r(z + y + e - 2)r(2 - e - z)r(2 - e - y) 

T(x)r(y)T(A-x-y-2e) ' 
r(2 - x - z - e)r(2 - j/ - z - e)r(-2 + x + y + ^ + e) 

r(x)r(y)r(4 - x - y - 1 - 2l) 



(2.15) 



(2.16) 

(2.17) 
(2.18) 



Going through the same step for the cut in Figure []J easily confirm the result we have 
derived here. 



2.2 Calculation with adjoint generators 

The same calculation outlined in the previous section can be described in a very efficient 
way using the representation of one-loop amplitudes in terms of adjoint generators [16] 
described in f ]A.3j) . We consider again a two-particle cut with a tree-level form factor and 
a one-loop four-point amplitude 

A w (l,2,h,l 2 ) = Yl TriF^-.-F^A^ax,...,^), (2.19) 

where (a\, . . . , cr 4 ) are permutations of (1, 2, Zi, Z 2 ). The sum contains three terms, as- 
sociated with the orderings (1, 2, Zi, Z 2 ), (l,2,Z 2 ,Zi) (1, Zi, 2, Z 2 ) . On the cut, we contract 
the colour labels a^, a/ 2 of the loop legs with a 5 a 'i a! 2 from the two-point tree- level form 
factor, as already described in the previous section. Doing this, the first two, and the last 
colour orderings will give rise to the following two colour traces 

fia h a h rj^./pai pas pa h pa h \ § a h a i 2 Tr(F 0l F° ! i F 0,2 F ai i ) , (2.20) 

respectively, where^l 

N 2 

Tr(F ai F a2 F a F a ) = AN 2 5 aia2 , £ a6 Tr(F ai F a F a2 F a ) = 4 5 aia-2 . (2.21) 

In the following formulae we actually divide these functions by a power of q 2 per loop. 
5 Because of the slightly unconventional normalisation of the adjoint generators, we have Ti(F a F b ) — 
2NS ab , and (F a F a )i, c = 2NS}, C . This different normalisation for the adjoint generators is responsible for 
the appearance of factors of 4 in (I2.20p . 
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Taking into account these relations, one instantly arrives at (12.61) . 

Before closing this section, we also quote two useful formulae when dealing with adjoint 
traces, namely 

Tr(F a AF a B) = Tr(F a A t F a B) + Tr(F a A)Tr(F a B) , (2.22) 
which, in the case where A 1 = —A, simplifies to 

Tr(F a AF a B) = -Tr{F a A)Ti:{F a B) if A 1 = -A . (2.23) 

This explains the factor of 1/2 in the second relation in (12. 2 1 p . 



3 Generalised unitarity for planar form factors 

In this section we compute higher-point form factors using generalised unitarity [16tfT7]. 
The explicit calculation will focus on the MHV three-point two-loop form factor, 

F 3 MHV ' (2) (1,2,3) := (0i 2 (pi)0i 2 (p 2 )^ + (j9 3 )|Tr(0 12 12 )(O)|O) , (3.1) 

for which we will now derive the complete result. Recall that the operator Tr(0i 2 0i 2 ) is a 
half-BPS operator and is the lowest component of the stress-tensor multiplet. Hence, we 
expect the result to be the same if this operator is replaced by any other operator in the 
stress-tensor multiplet [TJ. A different choice of operator usually requires a different choice 
of external states as well. We also notice that in TV" = 4 SYM there is nothing special 
about the particular choice of helicities of the external states in (13. ip - any other choice 
with the same total helicity and .R-charge are possible. In fact, we could have worked in 
a more invariant way by using supersymmetric form factors as introduced in [7J, and in 
the case of MHV form factors all helicity dependence factors out into the MHV tree-level 
(super) form factor, 

F 3 MHV ' (2) (1, 2, 3) = F 3 MHV,(0) (1, 2, 3) gf\l, 2, 3) , (3.2) 

where the helicity-blind function Q\ (1, 2, 3) depends only on the momenta through si 2 : = 
2pi • P2, s 23 ■= 2p 2 • p 3 , 8 3 i := 2p 3 ■ pi, and q 2 := s n + s 23 + s 31 ; finally, F 3 MHV ' (0) (1, 2, 3) 
is the tree-level (super) form factor. For the particular helicity configuration in (13 .ip this 
tree-level form factor is (12)/((23)(31)). It is easy to see that in all calculations this 
tree-level piece factors out naturally, and therefore the goal is to determine Q\ (1,2,3). 
We mention that this ratio function is also shared by the NMHV three-point form factor 
defined by factoring out the NMHV tree-level form factor, due to the chiral symmetry of 
three-point case [7J. We also note in passing that three-point form factors with adjoint 
fields always have a colour factor proportional to the structure constants f a ^ a ^. 

We begin by reconsidering the cut we are familiar with from the two-particle two-loop 
form factor, namely the two-particle cut in the q 2 channel, where the cut integrand is 
given as the product of a tree-level two-point form factor and a one-loop amplitude. As 
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we have seen in the two-point case, the double-trace non-planar amplitudes also contribute 
to planar form factors because of the colour index contraction 



Tr(T ai iT a; 2) [Tr(T a ^T a ' 2 ) Tr(T ai • • -T a ")] = N 2 Tr(T ai • • ■ T a ") . (3.3) 



Using the general formulae for the one-loop amplitudes (IA.lj) -f lA72~]) . we immediately ob- 



tain the corresponding cut integrand of the planar n-point form factor, 

F^(q 2 ) = [dUPS(h,l 2 ;q) UA n]1 (h,l 2 ,a(l),...a(n)) (3.4) 

q 2 — cut J L 

[n-lJ/2 

+ 2 ^ A^tih, ct(1), . . . , a(i), l 2 , a{i + 1), . . . , a(n)) 

i=l 

+ A n . x {l x , ct(1), . . . , a(n/2), / 2 , a(n/2 + 1), . . . , a(n)) 

where the last line contributes only when n is even. In particular, for three-point form 
factors we have 

ifV) 2 = fdLlPS(h,l 2 ;q) (3.5) 

q 2 — cut J 

[4A 5;1 (/ 1 ,/ 2 ,ct(1),ct(2),ct(3)) + 2A 5;1 (/ 1 ,ct(1),/ 2 ,ct(2),ct(3)) . 

Note that in the above presentation we have fixed the order of the external on-shell states, 
and stripped off the colour factor jV 2 Tr(T ai T a2 . . . T a "). We call this quantity a colour- 
stripped form factor, or colour-ordered form factor. The complete planar form factor can 
be expressed in terms of these colour-stripped form factor times the appropriate colour 
factor, and summing over non-cyclic permutation of external on-shell states. It is this 
colour-ordered form factor that we will focus on in the later calculation. 

We would also like to point out that colour-ordered form factors can also be computed 
directly using colour-stripped Feynman rules, which are simple generalisation of the well- 
known colour-stripped Feynman rules for amplitudes. The relevant diagrams are planar 
diagrams - in the sense that we fix the ordering of the external on-shell particles, while the 
off-shell momentum q of the gauge-invariant operator can appear in any possible position. 
Hence, in momentum space, diagrams can look non-planar, while in fact in colour-space 
they are planar. 

We will now focus on the computation of the three-point form factor. Unlike the 
two-point case, it is technically more difficult to lift the two-particle q 2 cut directly to the 
two-loop form factor, because of the presence of the one- loop five-point amplitude (13.51) . 
However, we have used these cuts and two-particle cuts in other channels as important 
cross-checks of our final result. 

For n-point form factors, not much is known about the set of two-loop integral func- 
tions appearing (except for n = 3 [H^ED]), and since dual conformal symmetry is not 
a feature of form factors it cannot be used to reduce the set of allowed integrals. Fur- 
thermore, planar form factors, as already seen in the case of the Sudakov form factor, 
contain non-planar integral topologies. Therefore, we follow a two-step procedure that 
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Figure 2: Double two-particle cuts and three-particle cuts of three-point form factor. 



does not rely on the knowledge of a particular basis of integral functions: in the first step 
we use double two-particle cuts to identify a basis. This set of cuts leads to sufficiently 
simple integrands so that their tensor reductions can be performed with relative ease and 
a simple set of integral functions appears. In this step we determine all possible integral 
topologies together with their coefficients, but are left with ambiguities in the numerators 
of the form ~ If, where U is one of the cut momenta. In the second step we perform 
more intricate and stringent triple cuts to fix any remaining ambiguities. One additional 
complication here is that sums over the helicities of the internal particles have to be per- 
formed. However, this can be dealt with elegantly using superamplitudes and super form 
factors [7J. This procedure leads to a unique answer that passes all consistency checks. 

More concretely, following the strategy outlined above we determined the integrals by 
considering the set of double two-particle cuts and three-particle cuts shown in Figure El 
These cuts are sufficient to determine the three-point form factor completely. We have also 
considered various two-particle cuts with one-loop form factor times tree-level amplitudes, 
and tree-level form factor times one-loop amplitudes which provided important additional 
cross checks, but did not lead to new integral functions. 

We now provide some more details on the procedure we followed. Starting from the 
two-particle cut expression f 13 ■ 5 j) . one can apply a further two-particle cut to the one-loop 
five-point amplitudes. The cut integrand is then given by the product of a two-point tree- 
level form factor and two tree amplitudes. We also consider the cuts which are given by 
a three-point tree form factor and two tree amplitudes. These types of cuts are depicted 
in the first two lines of Figure |2j The cut integrands are simple enough to perform the 
necessary tensor reduction directly. In this way, we find a set of simple integral functions 
with simple coefficients containing all integrals that appear in the final answer, given in 
Figure EE 

However, we are left with certain ambiguities due to l\ 2 terms in the numerator of 
integrals such as DBox and NBox in Figure [6j where Z 1)2 correspond to cut propagators. 

6 There are also additional non-planar integrals arising from the reduction of the double trace term 
in ()3.5|) . some of which are not allowed by the colour structure. However, these are all canceled after 
performing cyclic summation, and do not appear in the final result. 
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Figure 3: A particular three-particle cut of the three-point form factor. 



Such terms are not detected in the double two-particle cuts considered. Besides, there are 
also integrals which are not detected by these double two-particle cuts. Both problems 
can be fixed by considering three-particle cuts. 

The three-particle cuts on their own involve several integral topologies^ The cuts we 
have considered are shown in the third line of Figure |2j and involve up to six-point NMHV 
amplitudes and five-point NMHV form factors. The cut integrands are therefore much 
more complicated compared to double two-particle cuts, which makes it much harder to 
perform the tensor reduction directly in order to obtain a set of simple integral functions. 
That is why we chose to use first double two-particle cuts to write down an ansatz, which 
we then verify and refine using three-particle cuts. 

For these checks we do not need the reduction of the triple-cut expressions as we can 
make analytic comparisons of the integrands arising from the triple cut and those coming 
from the ansatz directly, by choosing a basis for spinors and expressing both integrands 
in this basis. 

To be more explicit, let us consider the particular three-particle cut shown in Figure [3] 
The cut integrand is given as a product of a four-point tree form factor, and a five-point 
tree amplitude. There are two kinds of contributions to the cut integrand, depending on 
which one is next-to-MHV: 

J dVdWSfo [^ 4 MHV,(0) Hi, -h, -h, 3) ^ MHV ' (0) (1, 2, h, k, h) 

+iF NMHV )( 0) Hi) _ K _ /3; 3) A MHV )( 0) (lj % ^ ^ ^ ^ 

where in the last equation we have used supersymmetric amplitudes and form factors [7] 
in order to perform the sum over internal helicities efficiently The fermionic integration 
can now be performed easily, and after switching back to component amplitudes and form 
factors the result is 



J -MHV,(0) (1;2)3) 



(12) (23) (31) 4 NMHV,(o) n+ 9+ ,_ 

immmw 5 >*3>*2>*ij 

(23)(31) p NMHV,(0),_ 7 - _,- o + \ 

r 4,SD I n > l 2 i l 3 i ° ) 



(2k)(kl 2 )(l 2 h){hl) 



(3.7) 



Note that the amplitudes and form factor in the bracket are the bosonic components with 
fixed helicities. i^^, HV '^ is the form factor with an insertion of the operator Tr(Fg D ) + . . . 



integrals which are simple products of one-loop integrals are not detected by three-particle cut, but 
they are ruled out by the double two-particle cuts. 
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Figure 4: There are twenty integrals contributing to the three-particle cut example dis- 
cussed in the text. 



(the dots refer to additional cubic and quartic terms) which is part of the stress-energy 
multiplet. The appearance of this form factor is due to the fact that both ^ r ^g D HV ^°' ) and 
j-.nmhv,(o) gkQj-g same NMHV factor, i.e. 

r NMHV,(0) r NMHV,(0) 

t-,NMHV,(0) _ "^SD _ ^4 /o o\ 

4 — t-MHV^O) ~~ r MHV,(0) ' V ) 

4,SD * / "4 

where the former is the super form factor with an insertion of the Lagrangian Ti(Fg D ) + - ■ ■ 
and the latter with Tr(0 2 ) \7\. By performing the fermionic integrations we extract the 
TJi 1 Vi 2 vt 3 component, and obtain F^^ V ^°\—l^,—l2,—l^,3 + ). This form factor can be 
calculated using BCFW recursion relations as in [7], with the result 

f nmuv,(o) ( +) _ (3|p 12 |4] 3 (l|p 23 |4] 3 

4 ' SD 1 ' ' ' } " S412 [41][12](3b 41 |2] + S2 34[23][34](l|p 3 4|2] 1 J 

(13) s 1234 



S34l(34)(41)(3|p 4 l|2](l|p34|2] 



We compared this cut integrand with the expression obtained from the three-particle cut 
of the integral presentation of the final result given in Figure El As shown in Figure HI 
there are in total 20 contributing cut integrals. Remarkably, this cut involves all integral 
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(a) 



Figure 5: The two appearances of the integral TriPent in the three-particle cut. They give 
different contributions to the cut integrand. 

topologies appearing in Figure [6] and we find complete agreement with the above non- 
trivial cut integrand_|. This provides one of the strongest consistency checks that our 
result summarised in Figure E] is indeed correct. 

There is an important subtlety we would like to point out. As we mentioned before, 
we used double two-particle cuts on (13.51) in order to first detect all the integrals which can 
appear. We find that TriPent and DTri actually receive two identical contributions from 
the single-trace and the double-trace term in (13.51) . We have drawn the two contributions 
from the planar and non-planar part in Figures [5](a) and \Mb), respectively. Similar com- 
ments apply to DTri. This observation is relevant when we consider unitarity cuts. For 
example, in the three-particle example we considered above, two different presentation of 
TriPent give a different contribution to the cut integrand - indeed, as shown in Figure 
El the positions of h and h are swapped, and therefore the total contribution to the cut 
integrand is 

1 9 2 Sl2S23 1 Q 2 Sl2S2Z 

2 (h + P2) 2 (h - P3) 2 (h + k) 2 (h + q) 2 + 2 (h + P2) 2 (h ' P3) 2 (h + k) 2 (k + q) 2 ' ( ' } 



We conclude this section by presenting our complete result for the two-loop three- 
particle form factor in M = 4 SYM obtained from generalised unitarity cuts: 

2 

qP = ^(DTrii + DBoXi) + TriPent + NBox + NTri + cyclic , (3.11) 

8=1 

where in Figure [6] we have presented compactly all integrals, including their precise coef- 
ficients and numerators some of which are loop dependent. 

The integral functions appearing in (13.111) have six or seven propagators and notably 
both planar and non-planar topologies appear, some with irreducible, loop-dependent 
numerators. The maximum power of loop momentum appearing in the numerators is 
one. We will now proceed to discuss the numerical evaluation of these integrals. 

8 Note that no parity-odd terms can be present in the final result for the three-point form factor since 
there are not enough independent momenta to form a non- vanishing contraction with e^ vp \. For the 
same reason we also expect there should be no /x 2 -terms in the expressions of the cuts, and hence it is 
sufficient to use four-dimensional unitarity. 
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3.1 Evaluation of the integral functions 



Several of the integral functions appearing in our result (13. lip have been computed and 
explicit expressions can be found in [T8H20] - in particular all the six-propagator integrals 
are known. For the remaining integrals TriPent, DBoXi and NBox explicit analytic 
expressions are not known, in part due to the peculiar numerators we found. In principle, 
there exist algorithms that allow to reduce these integrals to known master integrals 
[T8H20] . However, in this work we chose to use an independent route, and expressed them 
in terms of Mellin-Barnes (MB) representations (see [21] for a pedagogical introduction) 
and used well-established codes [22] to evaluate them numerically to high precision. 

The MB form of the planar integrals can readily be constructed with public pro- 
grammes like AMBRE [48, 49J and we will not quote their explicit forms here. However, 
the MB representation of the NBox integral cannot be found with AMBRE due to well- 
known issues of this programme with non-planar topologies with several scales. Therefore, 
we constructed an MB representation of NBox directly from its Feynman parameter form. 
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The result is an eight-fold MB representation of the form 



2(27ri)*r(-l - 3e) J l\ {dZlL { Zl)) U V W X 

T(-e - z u )T(-e + z 4 )r{l + Zi 3456 )r(l + z 157 )T(-l -e + z 3 -z s )x 

r(-2 - 2e - z x - ^6s)r(-2 - 2e - 2i 34 57 8 )r(-2 - 2e - 2123467s) x (3.12) 

r(-2e - 23 + ^s)r(l + gi68)£(l + 2 2 7s)r(3 + 2e + zi 2345678 ) 

r(-2e - z 3 )r(-l -2e-Z3- z^)V(-l - 2e + z 34 - z B )T(-2e - 234 + * 8 ) ' 

where we have introduced the shorthand notation — z % + z j + • • • + z k-> an d 

s 12 S 2 3 S31 

m =: — , u := — , w := — . (3.13) 
q 2 q 2 q 2 

Note that for sake of brevity we have dropped here the terms of the numerator which 
are linear in loop momentum t\ they lead to a number of similar eight-fold MB integrals. 
Furthermore, due to the T(— 1 — 3e) denominator the integral effectively becomes seven- 
dimensional |2T| . In this sense this integral is the most complicated and numerically the 
most challenging contribution to the form factor, since the planar topologies DBoXi and 
TriPent require at most three- and four-fold MB integrals, respectively. 

We have evaluated (13. lip by expressing all six-propagator integrals by their known 
analytic formulae and using MB . m by [22] for numerical evaluations of the MB representa- 
tions of all seven-propagator topologies. The result is a power series starting as e -4 which 
we have computed up to and including finite O(e ) terms. We present here a few results 
of our numerical evaluation at four kinematic points (— S12, — s 23 , — s 31 ): 

(1.1.1) : ^ + | + ^-^?- 18.2484 ± 0.02 + O(£ ), (344) 

(1 . 1. 2) : k ~ ^ + ™ " — " T-3182 ± 0.02 + 0( £ ) , 

e 4 e 4 e 2 e 

(1 . 2. 2) : « - + - ^ + L868 6 ± 0.02 + 0(e) , 

(1 . 2 . 3 ) : k ~ ^ + ^ " ^ + I" ± 0.03 + O(e) . 

e 4 e A e 2 e 

Here we have only quoted the errors of the finite terms. The numerical error of the 1 /e 4 
and 1/e 3 are negligible, while for all kinematic points we have investigated, a conservative 
estimate of the errors of the 1/e 2 and 1/e terms is 10 -13 and 10~ 7 , respectively. 

In the following we will study this two-loop result in detail and, motivated by the fact 
that form factors have universal collinear limits, we will define a finite remainder very 
much in the spirit of [UI2]- We will then determine the symbol [23] of the remainder 
function, and from that derive its analytic expression, which we will compare to our 
numerical results. 
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4 Exponentiation of the form factor and the remain- 
der function 



In this section we want to consider the possibility that higher-loop form factors in M = 4 
SYM obey a similar exponentiation relation as MHV loop amplitudes [UI2] for small 
numbers of on-shell particles, and understand from which number of particles a remainder 
function has to be added. We expect this remainder to appear for a smaller number of 
particles compared to amplitudes because only Lorentz symmetry and dilatations are 
unbroken (up to infrared divergences). This leaves us with 3n — 7 parameters on which 
a remainder might depend. Indeed we will find that the three-point remainder function 
depends on two variables. Similar observations were made at strong coupling for (1 + 1)- 
dimensional kinematics in [5], where a non-trivial remainder appears at four points. 

For the Sudakov form factor, exponentiation is evidently true even in QCD, and in 
J\f = 4 this was explicitly proved in [8]. Specifically, one finds 



F^\ q \e)-\{F^\ q \e)) 2 = (-g 2 )^ 



% + ^ + 0(e) 
e z e 



(4.1) 



where is the L-loop Sudakov form factor. The result in ( 14. ip was re-derived in 
appendix B of [6], where it was recast in a slightly more modern language as 

F {2 \q\e)-\(F^\ q \e)) 2 = /< 2 > (e)F« (g 2 , 2c) + C< 2 > + 0(e) , (4.2) 

where /< 2 )(e) = / (2) + f {2) t + /fe 2 , with 

fP = -2( 2 , /« = -2C 3 , (4.3) 
and with one relation between and C^- 2 \ namely^ 

C <2) = h_ + l_ (4 _ 4) 

2 18 

For comparison with amplitudes, we recall that the four-point MHV amplitude (divided 
by the tree-level amplitude) satisfies [U|2] 

M (2) (e) - liM^ie)) 2 = m p (e)Ml\2e) + C% + 0(e) , (4.5) 
with / £Ue) = /o (2) + e/f ) + e 2 /f,and3 

,(2) f (2) 

7o (2) = -C 2 = , A (2) = -Cs = 4r • (4-6) 



9 We note that / 2 (2) will be fixed later in ([420]) fr om collinear factorisation. 

10 As observed in [6], the factor of 1/2 in the result (|4.6|) is a matter of convention - it can be understood 

(2) (2) 

once one recalls that /q and are written in a convention where the t Hooft coupling obds is twice 
as that used in [SUB]. It is the combination a/( 2 ' which must be independent of any conventions used to 
define the coupling. 
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4.1 Collinear factorisation 



In this section we would like to discuss collinear factorisation and the exponentiation of 
infrared divergences of form factors. This will lead us to the definition of a finite, scaling- 
invariant and regulator-independent remainder function with trivial collinear limits. 

We begin by recalling that an important hint that MHV amplitudes in M = 4 may 
exponentiate came from the study of collinear limits [1] where two adjacent momenta 
a and b become parallel. It is well known that scattering amplitudes have a universal 
collinear factorisation behaviour, which is governed by splitting amplitudes [T4"lll~5ll5T)] . 
These quantities only depend on the helicities and the momentum fractions carried by 
the two legs becoming collinear, but are completely blind to other details of the process. 
If we consider M = 4 SYM and focus on the case of MHV amplitudes, the helicity-blind 
ratio /An '' = Mn(pi, ■ ■ ■ ,p n ) therefore obeys the following factorisation: 

L 

M n L) (p u ...p a ,p b ,... ,Pn)^J2 Mi ^ Pl > ■■■,Pa + Pb,-- .,p n )r (L - l) (e; z, s ab ) , (4.7) 

1=0 

where = 1, M.n^ = 1, and 



1 - 2 F X I 1, -e, 1 - e, — 1 - 2 F X I 1, -e, 1 - e, 



(l)f \ C V ( ~ s ab 

T y€, Z, S ab ) . I - j i i - ' i i 

€ z V J \ Z J \ Z — 1 

(4.8) 

to all orders in the dimensional regularisation parameter e [511 152] 1^1 

For general large- N gauge theories, a unitarity-based proof of (14. 7p at any loop order 
has been given in jl5] . The important point we would like to make is that the main steps 
of that proof apply directly to form factors as well, which therefore share the universal 
factorisation properties of amplitudes, including (14. 7p . which we write as 

L 

G n L) (Pl, ■ --PaiPb, ■ ■ ■ ,Pn)— * ■■■iPa+Pb,-- ■ , Pn)r {L ~ l) (e; Z, S ab ) , (4.9) 

1=0 

where we have defined 

Qi L) ■= W° • (4-10) 

We stress that the splitting amplitudes appearing in (14. 9 p are the same as for amplitudes^ 

Remarkably, loop splitting amplitudes in M = 4 SYM obey a cross-order relation 
which at two loops is [1] 

r (2) (e; *, s) = ~ (r«(e; z, s)) 2 + f 2 \e)r i - 1 \2e ] z, s) + 0(e) , (4.11) 



n The all-orders in e expression for r^'(e; z, s a b) in (|4.8j) was also rederived in [53] using one-loop MHV 
diagrams. 

12 We have checked this for bilinear operators following the proof of [5S], however we caution the reader 
that the proof presented in that paper does not extend automatically to operators such as Trtp n for n > 2, 
and collinear factorisation should be re-examined carefully in these cases. 
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which has generalisations to higher loops [2]. Importantly, this implies that the correct 

f2l 

collinear factorisation of the two- loop amplitude Ai n is captured by the combination 

M<$ := \{M^\e)f + f 2 \e)M^{2e) , (4.12) 
which is given purely in terms of the one-loop ratio functions and f^ 2 '(e). The only allowed 

- — (2) (2) 

discrepancy between Ai n and the full two-loop ratio Ai n would be a function that is 
finite in all collinear limits. As is well known, in M = 4 SYM this difference turns out to 
be a constant for n = 4, 5 and a non-trivial function n > 5 - the remainder function. This 
behaviour is linked to the fact that loop amplitudes in M = 4 SYM obey an anomalous 
dual conformal Ward identity [54T. I55] , of which A4n ^ is a particular solution. For n = 4, 5 
no dual conformal cross-ratios exists, and the only homogeneous solution is a constant, 
while for n > 5 there are 3n — 15 cross ratios and explicit calculations have confirmed the 
necessity to add a remainder function TZ n (ui, . . . , W3 n _i 5 ) 



Given the similarities between amplitudes and form factors, specifically concerning 
their collinear factorisation properties and the universality and exponentiation of their 
infrared divergences [361I56T459] , it is natural to expect that loop form factors should have 
a similar exponentiated form. Let us investigate this point in more detail. 

To begin with, we have checked explicitly in a specific example our claim based on 
a simple generalisation of [45] that form factors enjoy collinear factorisation properties 
identical to those of amplitudes. Specifically, we have confirmed that the one-loop MHV 
form factors derived in [6] indeed obey the anticipated collinear factorisation (14.91) . For 
concreteness we now consider the one-loop three-point form factor ratio function (see 
(3.17) in ©), 



su)- 6 + (-s 23 r e + (s 31 )- e (4.13) 
+{-S23)-h( - ^) + (s 31 )-*h( - *>) - (-q 2 r £ h(- 

\ Sai/ V Sr>?J 



s 31 7 v s 23 7 S23S31 

\ si 2 / V s 31 J V S31S12/J 

with h(x) = 2-^1(1, — e, l — e,x) — l and q 2 = (pi + p 2 + P3) 2 = S12 + S23 + «3i- 

In the collinear limit 1||2 we have S12 — > 0, q 2 — > S23 + S31. Setting, in the limit, 
s 3i/q 2 z an d S23/V — 1 — 2, we find 



g 8 W rM(e;s l2 ,z) - ^(-(P + p 3 ) 2 )^ , (4.14) 



with P — pi + pi- Note that the splitting amplitude comes from terms 1, 4 and 11 in 
(I4.13p . and the second term is nothing but the one- loop two-point (Sudakov) form factor. 
Reproducing the expected factorisation requires a delicate conspiracy among various terms 
in ([mi). 
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4.2 Iterative structure at higher loops 

Having established the correct collinear factorisation at one loop, we can now write down 
a general two-loop MHV form factor ansatz as follows, 

Si 2) = + / (2) (e) <#>(2e) + , (4.15) 

where the ratio is defined in (14.10p . Note that we allow for a potential form factor 
remainder function lZ n on the right-hand side of (14. 15ft . A few key properties of the form 
factor remainder are: 

1. It must free of infrared divergences, and 

2. It must be finite in all collinear limits. Upon properly normalising the remain- 
der function by adding to it an n-independent, transcendentality-four constant, we 
expect that the normalised remainder behaves as 

K n -> ft n _i , (4.16) 

in a simple collinear limit. Furthermore, we also expect that 

3. It must be expressed in terms of transcendentality-four functions, of which we will 
shortly discuss the symbol. 

4. It is a rescaling invariant function and hence it depends on Mandelstam variables 
only through their ratios. But due to the lack of dual conformal (super) symmetry 
of form factors it does not enjoy the more restricted dependence on conformal cross 
ratios of amplitude remainders in planar M = 4 SYM. For this reason, a non-trivial 
form factor remainder can already appear at n = 3. 

The first property follows from the exponentiation of infrared divergences. Pleasingly, 
we have checked it explicitly in the case of the three-point, two-loop form factor. 

The second property follows from the analysis of collinear limits at higher loops |28j . 
This analysis was performed in the case of amplitudes in Section 7 of [60], following 
[28] . and it is instantly extended to form factors because their collinear factorisation, as 
expressed in ( 14. 9p . works in the same way as in the amplitude case ( 14. 7p . In particular, 
using ( 14. 9p . one finds that under a simple collinear limit, the scalar functions Git and 
Qn^ behave as 

G$> 0£i + r«(e;* jaa6 ) , (4.17) 

-> + ^(e; *, + r {2) (e- z, s ab ) . 

Using the fact that splitting amplitudes obey an iterative formula identical to the homoge- 
neous form of the BDS conjecture for the amplitude [T], expressed by (14. lip , we conclude 
that under a simple collinear limit, 

G%Xe)-\(G£\e)) 2 -fV\e)GV{2e) 

-> " - / (2) (e)^ 1 (2e) . (4.18) 
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Equation (I4.15P defines the form factor remainder TZ n as 

^ ■= G$>1?) " l(SP{e))'-f®{e)gV(2e) - + O(e) , (4.19) 

and it follows from (I4.18P that in the simple collinear limit TZ n — > TZ n -i, as anticipated 
in f)4.16p . The outcome of this brief analysis is that writing the L-loop form factor 
as an L-loop BDS-like ansatz plus an L-loop remainder gives a very efficient way to 
organise the presentation of the result - the BDS part captures infrared divergences and 
simple collinear factorisation, while the remainder, is finite in four dimensions and in the 
collinear limit flows smoothly from the n-point to the n— 1-point remainder. An important 
comment is in order here. The iterative structure in (I4.18P is driven by the iteration of the 
splitting amplitude ( 14. lip . For this reason, the function /™(e) is completely determined 
to be the same as in the amplitude case. In our normalisations for the coupling constant, 
we have 

/< 2 >(e) = -2( 2 - 2C 3 e - 2( 4 e 2 . (4.20) 
In turn, this implies that we can determine the constant introduced in the iteration 

_____ tiy\ ______ 

of the Sudakov form factor, see (jOl . Plugging f^> = -2( 4 into (jOjl . we find 

C (2) = 4( 4 • (4.21) 

Importantly, the constant is independent of the number of legs. The consequence 
of this observation is that the three-point form factor remainder function as defined in 
(I4.19P (for n = 3) with the value for just found, must have a very precise collinear 
limit: 

nf ] ->• , (4.22) 

in any of the three possible simple collinear limits. The third property is simply a conse- 
quence of maximal transcendentality. 



4.3 Numerical results for the three-point remainder at two loops 

As explained in the previous section, it is natural to define a remainder function, which 
at three points is given by 

Tlf := gf\e) - \{^\e)f - /«(e)C#>(2 e ) - C® + 0(e) . (4.23) 
This quantity is finite, regulator-independent and vanishes in all soft and collinear limits. 

(2) 

Our numerical evaluations confirmed that all infrared- divergent e~ p terms in 7^.3 do 
indeed cancel: for p = 3, 4 numerical errors are negligible, while for p = 2 and p = 1 we 
have confirmed this up to about 10~ 13 and 10~ 7 , respectively. This is a stringent check 
that our result ( 13. lip obtained using generalised unitarity is indeed correct. We have 
collected in Table 1 a few values of the remainder function obtained with our numerical 
programmes. The remainder is defined in (I4.23P where we have set C^ 2 ' = 4£ 4 and 
fp = -2( 4 . We will later compare this numerical remainder to our analytical expression, 
see Table 2. 
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(u,v,w) 


numerical l-C^ 


est. error 


(1/3, 1/3, 1/3) 


-0.1519 


0.02 


/ 1 / A -1 1 A -1 /C~\\ 

(1/4, 1/4, 1/2) 


-0.1203 


0.02 


(1/5, 2/5, 2/5) 


-0.1301 


0.02 


(1/2, 1/3, 1/6) 


-0.1080 


0.03 



Table 1: Some numerical values of the remainder Tlf defined in (14.231) . with = 4^ 
and /( 2 ) = — 2£ 2 — 2( 3 e — 2£ 4 e 2 . Infrared poles cancel with negligible numerical errors. 



4.4 The symbol of the three-point remainder at two loops 

The question we would like to discuss now is how we can constrain the form factor 
remainder function nffl defined in (14.191) using its symbol [23] @ We will focus on the 
case of the three-point form factor calculated explicitly in this paper. The conclusion of 
the analysis presented in this section is that, in the three-point case, there is a unique 
non-trivial solution for the symbol of the remainder function (up to an overall constant) 
to be presented shortly. 

In order to write down an ansatz for the symbol of the remainder function, we make 
the following physical considerations which lead to important constraints on its structure: 

1. We should understand which variables the remainder function can depend upon. 

2. The symbol must satisfy a certain first-entry condition [30], arising from the require- 
ment that the two-loop form factor remainder must have discontinuities only across 
the physical branch cuts starting at Pj = 0, where Pj's are appropriate kinematic 
invariants made of sums of external momenta (in the three-point case considered 
here, these could be Si 2 , s 2 3 and Si 3 ). This implies that the symbol of the 
two-loop remainder must be of the form 

5(2) = Pj ® S\vf] , (4.24) 

disc j 

where is the discontinuity of the two-loop remainder in the Pj-channel, and 
the sum is extended to all channels where there is a physical discontinuity. 

3. The second and the last entries of the symbol must also satisfy certain additional 
constraints, to be discussed below. 

4. The symbol must have trivial collinear limits. 

5. The symbol must satisfy all relevant symmetries. 

6. The ansatz must be the symbol of a local function, which in turn means that the 
symbol must satisfy certain integrability conditions. 

13 Similar discussions have appeared for the symbol of amplitude remainder functions in [30H35) . 
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We now focus on the three-point case and discuss the implementation of the properties 
listed above. 

1. To begin with, we wish to claim that the symbol of the three-point form factor 
remainder function takes its entries only from the following set of six kinematic variables 

{u, v, w; 1 — u, 1 — v, 1 — w} , (4-25) 

where u, v and w are the scale invariant ratios defined in (13.131) . This statement is highly 
non-trivial and is an important constraint, as it does not allow for dependence on any other 
variables such as 1 + u, or 1 + uv etc. Our claim is substantiated by explicit inspection of 
the two-loop master integrals with one off-shell and three on-shell legs [T9l[20] in terms of 
which our two-loop, three-point form factors is expanded. We have calculated the symbol 
of the relevant integral functions, with the result that they always depend on the variables 
listed in ( 14 . 2 5 f) only and that the first entry of the symbol is always taken from the subset 
{u,v,w}. 

2. The first-entry condition requires that the first entry of the symbol must be either 
u, v, or w, but the remaining three variables 1 — u, l — v, l — w cannot appear as possible 
first entries. This is also confirmed by inspection of the integral functions appearing 
in [19l[20] as mentioned under 1. 

3. Furthermore, the second entry of the symbol is restricted in the following way |30j . 
If the first entry is u, then the second entry can be 1 — u or v or w, but not u or 
1 — v or 1 — w; similar conditions in the case when first entry is v or w are obtained 
by cyclic permutations of u, v and w. Finally, we impose a last entry condition on the 
symbol [3"4l3~T] . by requiring it to be always of the form uj (1— u), or vj (l—v), or w /{l—w). 
If we insist in using the variables in (14.251) as entries for our symbol, the above conditions 
relate pairs of coefficients, namely the coefficient of a term in the symbol whose last entry 
is u is the opposite of the coefficient of the term whose first three entries are the same, 
with the last one equal to 1—u. And similarly for the pairs v and l — v, and w and l — w. 

At this point, we observe that imposing the first three conditions, one obtains an 
ansatz for the symbol in terms of 324 independent parameters. 

4. Next, we recall that the three-point remainder must vanish in all collinear limitsl^l 
Imposing this constraint term by term on the symbol is too restrictive, as there can be 
cancellations among different terms. Therefore, we keep the most general ansatz and 
impose that it vanishes in any of the three simple collinear limits. We also recall that, in 
order for a symbol to vanish, it is enough that one of its entries smoothly goes to one in 
the limit. 

For definiteness, let us focus on the s i2 — > 0, or u — > limit, which we take as 
follows: we leave untouched every entry equal to u in the symbol; however we can replace 
1—u 1, and furthermore, since v+w = 1 in the collinear limit, we can replace w — > l—v. 
Requiring the symbol to vanish in this limit, we obtain three sets of equations, for the 
three independent simple collinear limits s\2 — > (or u — > 0), S23 — >■ (or v — > 0) and 
S31 —^0 (w — > 0). Altogether, the collinear equations impose 99 independent constraints 
on the symbol. 

14 We would like to thank Paul Heslop for a very useful conversation on collinear limits of symbols. 
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(2) 

5. Next, we impose that the remainder TZ 3 must be totally symmetric in u, v, w. 
This leads to 52 linear relations between the coefficients. 

6. Finally, we impose the integrability condition on the ansatz, namely the constraint 
that the symbol is derived from a local function. This is implemented as follows [SUED]: 
one takes the symbol 

S = J^^i ® ■ ■ • <g> w n , (4.26) 

and replaces two consecutive entries Wi, Wi + i by dwi A dw. l+ i. The resulting expression 
must vanish: 

dwi A dw i+ i (wi ® ■ ■ • ® Wi-i ® u7j + 2 <8> • • • ® = . (4.27) 

There are three different integrability conditions we have to impose in our case0 Re- 
markably, by requiring these to be obeyed, we find that there is a unique solution for the 
symbol of the two- loop remainder (up to an overall normalisation). It is given by the 
following compact expression: 

S {2) = -2u®(l-u)®(l-u)®- — -+u®(l-u)®u®- — - 

u u 

/-, x 1—v /A . 1 — w 
— u ® (1 — u) ® v ® u ® (1 — w) ® u> ® 

/-I \ 1 — v /-, \ 1 — It 

— M ® V ® (1 — U) ® « ® 1> ® (1 — V) ® 

V u 
1 — u 1 — v 

+U ® f ® W ® htl®t)®W® 

l — W .„ v l 1i) 

+u ®v®w® u ® w ® (1 — u) ® 

1 — "U I — f 

+M ® W ® U ® h M ® W ® V ® 

1 — W 1 — M 

+m ® u> ® ?; ® « ® it; ® (1 — w) ® 

+ cyclic permutations . (4.28) 



The next challenge is twofold: firstly, we wish to determine the function whose symbol 
is given by (I4.28P ; and secondly, we wish to determine terms missed by the symbol, 
e.g. terms of the form 7r 2 x F2 where F2 is a sum of transcendentality-two functions with 
rational coefficients. 

In this respect, there is an additional piece of information about ( 14 . 2 8 [) that we would 
like to mention. Our symbol defined in (14.281) satisfies an important symmetry 
constraint [61] discussed in [23], namely 

SfL - 41 - SSL + 41 - (« » c , d) = . (4.29) 

According to a conjecture of Goncharov, symbols with this peculiar property can always 
be obtained from a function involving logarithms and classical polylogarithms Li^'s with 
k < 4 only [6TJ[23]. The explicit solution we will present in the next section will confirm 

15 In general, there are p — 1 integrability conditions for a transcendentality p symbol. 
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this expectation beautifully. As we will show in the final part of this paper, there is an 
alternative way to obtain an analytic result of the form factor remainder in terms of two- 
dimensional harmonic polylogarithms [62J. This is due to a remarkable relation between 
the Af = 4 form factor and the planar, maximally transcendental part of the two-loop 
QCD amplitude for H — >■ ggg recently obtained in [i~2|[T3]. 



4.5 The analytic remainder function 

The remaining task now is to find a transcendentality-four function whose symbol is given 
by (I4.28p . Recall that the symbol only takes entries from the list {u, v, w, 1—u, 1—v, 1— w} 
and has the symmetry (I4.29P , which implies the result should be expressed purely in terms 
of classical polylogarithms of degree up to four and logarithms [5TJI23] • This however does 
not fix a priori the allowed arguments of these functions, but the arguments of individual 
functions must be such that the symbol of that function has only entries from that list. 
Taking these considerations into account, the most general ansatz will be built from the 
following set of functions: 

log x\ log x 2 log x 3 log x 4 , Li 2 (xi) logx 2 logX3 , Li 2 (xi)Li 2 (x 2 ) , Li 3 (xi) logx 2 and Li 4 (xj) , 

(4.30) 

where we found it sufficient to take the possible arguments X; from the list 



u, v, w, 1 — u, 1 — v, 1 — w, 1 



1 1 
- 1-- 1 

U V 



1 

w 



uv 

w ' 



vw 
u 



wu 

V 



(4.31) 



Imposing the constraint that the ansatz has the same symbol as (14.281) one can easily 
find a solution. We have then applied various polylogarithm identities to simplify the raw 
solution obtained in this way. The final result takes the remarkably simple and compact 
form 
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where U\ — u, it 2 = v and us 
J 4 (z) := Li 4 (z 



w and we have introduced the function 

2/ .A i„„3 



log(-, ) Li 3 (, ) + i^Li 2 (, ) -^ 
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(4.32) 



(4.33) 



It is curious to note here that J 4 (z) is almost identical to the function D 4 (z) introduced 
by Ramakrishnan. The functions D m (z), m > 2, are generalisations of the Bloch-Wigner 
functions (see [65] for an inspirational exposition of these topics and references). 

In the representation obtained above we have already taken into account beyond-the- 
symbol ambiguities which arise due to the fact that the symbol is blind to transcendentality- 
four terms of the form 7r 4 or n 2 x {logXjlogXj, Li 2 (xj)}. It is a simple task to fix these 
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ambiguities using constraints from permutation symmetry and collinear limits. In our case 
it was sufficient to add the £4 term to get a symmetric function, that is smooth throughout 
the Euclidean region defined asO<M<f,0<t'<f,0<w<f and u + v + w = I, and 
vanishes in all collinear and soft limits. 

Finally, we have collected in Table 2 results from our numerical evaluations in Section 
13.11 and compared them with the exact result (I4.32p . This also serves as confirmation of 
the overall normalisation of the remainder, which is not fixed by the symbol alone. 



(u,v,w) 


numerical TZ^ 


est. error 


analytic TZ^ 


(1/3, 1/3, 1/3) 


-0.1519 


0.02 


-0.148966 


(1/4, 1/4, 1/2) 


-0.1203 


0.02 


-0.134873 


(1/5, 2/5, 2/5) 


-0.1301 


0.02 


-0.136454 


(1/2, 1/3, 1/6) 


-0.1080 


0.03 


-0.125366 



Table 2: Comparison of numerical and analytic values of the remainder 71% 



4.6 Connections between amplitude and form factor remainders 

We have seen earlier that the three-point form factor remainder depends on three variables 
u, v and w defined in H3 . 131) . where u + v + w — 1. In (14. 28f) we have presented the symbol 
of this remainder. An a priori entirely different calculation is that of the six-point MHV 
amplitude remainder, whose symbol was explicitly evaluated in [23] using the results 
of [21]. This symbol can be expressed in terms of the three independent cross-ratios u, 
v and w one can write down with six lightlike momenta satisfying Yli=iPi = 0' as wen 
as three additional variables y u , y v , y w . An explicit compact expression for this six-point 
amplitude remainder symbol is given in Eq. (21) of [32] H It can be seen from that formula 
that the symbol is naturally decomposed into two terms, 

5 6 2 L P i = 4?L P iK^ w ) + si 2 L P i( u ' v ' w ->yu,yv,y w ) ■ (4.34) 

The second term in ( I4.34[) is given by [32J 

^6 2 L P i( M > v , w; y u , y v , y w ) = ~ [u sg) (1 - u) ® y u y v y w - 2u <g> v <g> y w ] <g> y u y v y w 

+ permutations. (4.35) 

On the other hand, it is straightforward to check that Sq amp i(u, v, w) is in fact identical 
(up to an overall multiplicative constant) to the symbol of our two-loop three-point form 
factor remainder (14.281) upon identifying the three (unconstrained) cross-ratios of the 
amplitude calculation with our three (constrained) variables u, v and w defined in (I3.13p . 

Note that 5g 2 ^ mpl (-u, v , w) has trivial collinear limits both when it is part of a form 
factor symbol, where u + v + w = 1, and when it is part of the six-point remainder 



16 The definitions of the y variables can also be found in Eq. (15) of [35]. We will not need their explicit 
expressions here. 
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symbol. In the latter case the three cross-ratios u, v and w are in general unconstrained, 
but in the collinear limit they become dependent - for instance, in the collinear limit 
w — > one also has that u + v — > 1 . 

This is reminiscent of a similar interesting coincidence observed at strong coupling 
in [5]. There, the four-point form factor remainder was evaluated in the particular case 
of (1 + l)-dimensional kinematics, and found to be expressible in terms of the eight-point 
MHV amplitude remainder in (1 + l)-kinematics, which was determined in [63, 64J. We 
note that the relations between form factor and amplitude remainders at strong and weak 
coupling seen in [5] and in the present paper can by no means be anticipated from the 
explicit calculation. For instance, in our weak coupling calculation we find that the three- 
point form factor remainder contains several integral functions - in particular non-planar 
ones - which do not appear in the corresponding six-point amplitude calculation. 

Incidentally, we note that while the four-point form factor remainder and the eight- 
point MHV amplitude remainder in of (l + l)-dimensional kinematics depend both on two 
independent variables^ in our case the number of independent variables is different - the 
form factor remainder depends on two independent simple ratios, and the amplitude's on 
three cross-ratios. 

Perhaps this connection between amplitude and form factor remainders is another 
indication that the paradigm "loop observable = rational coefficient x integral function" 
may obscure the simplicity of the final result of a loop calculation. Furthermore it is very 
likely that the symbol technology will allow to simplify the expressions of the (higher- 
loop) integral functions themselves, such as those appearing in (14.281) and in QCD. For a 
recent example in QCD see [55] . 

4.7 A surprising relation with QCD 

In this final section we wish to discuss an intriguing connection of our result with the 
recent work of [12]. In that paper, the two-loop helicity amplitudes for H — > ggg and 
H — > qqg were computed in the large top mass limit. In this approximation the top quark 
can be integrated out at one loop and produces a new effective vertex of the form Hgg. 
If we consider the process H — > ggg, which is the case of interest for our discussion here, 
the calculation is equivalent to that of a three-gluon form-factor with a TrF 2 operator 
insertion, where is the gluon field strength. As we noted earlier, the operator Tr0f 2 
- whose form factor we consider in this paper - is the lowest component of the stress- 
tensor multiplet, which among many other operators contains the on-shell Lagrangian 
£ . s . = Tr-Ff D + ■ ■ • , where Fg D is the self-dual part of the field strength, and the dots 
stand for cubic and quartic terms in the Lagrangian. Up to a universal helicity-dependent 
prefactor the form factor of the operator with three-gluons (g~ (pi)g~ (p2)g + (p3)\£o.sXQ)\ty 
is equal to the form factor we have been discussing in this paper, and on the other hand 
is related to the object considered in [T2] in QCD. 

The full QCD result of [12J is of course very complicated and has no resemblance to 

17 Recall that for this restricted kinematics, the 2n-point form factor remainder depends on 2n — 2 
variables, whereas the amplitude remainder depends on 2n — 6 cross-ratios. 
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(I4.32p . The story becomes more interesting if we focus on the planar part of the finite two- 
loop remainder A n defined in (5.8) and (5.20) of [12], and with explicit formulae in terms 
of 2dHPL's given in Eons. (B.l) and (B.7) of the same reference for the two different gluon 
helicity configurations.]^ It is important to note that these finite remainders are defined 
using the formalism of [36] widely used in QCD to remove the universal two-loop infrared 
divergences, and are different from the BDS finite remainders used in this paper [H[2]; in 
particular the former does not vanish in collinear limits. However, it is easy to relate the 
two definitions in J\f = 4 (see pQ), because one simply has to subtract ^(F^) 2 — 2( 2 F^ 
from Catani's remainder EB1 in order to arrive at the BDS remainder. 



The interesting observation arises as follows: First, take the finite two-loop remainders 
Aa or Ap from Eqns. (B.l) and (B.7) of [12], corresponding to the helicity configura- 
tions {g + g + g + ) and {g + g + g~) for the external gluons, respectively, and extract the terms 
containing functions of maximal transcendentality degree (four in this case). Then map 
them to BDS remainders using the prescription outlined above with replaced by the 
degree- two parts of the one-loop amplitudes Aa and AP (Eqns. (A.l) and (A. 4) in [12]), 
which coincides up to a factor of two with the e° term of our one-loop form factor in 
(14.131) . Finally, calculate the symbol of the resulting expression. Surprisingly, one finds 
exactly (14.281) . for both helicity configurations!^! This implies that our form factor re- 
mainder (14.321) is contained in the full QCD result of [12], and allows to present the latter 
in a more compact form. It would be interesting to see if symbols can be used effectively 
to simplify the terms of lower degree (transcendentality) in the expressions of [12], or in 
other QCD results. 

Alternatively, this connection provides us with an additional strategy to reconstruct a 
function from the symbol ( 14 . 3 2 [) by simply projecting out the functions of maximal degree 
and performing the subtraction described above to get the BDS remainder. The result 
obtained in this way is given in terms of 2dHPL's and the remaining beyond-the-symbol 
ambiguities can be fixed by considering symmetries and collinear limits. In this case the 
necessary correction terms can also be read off directly from e.g. A? in (B.6) of [T2J: one 
has to keep all terms of the form it 2 x F 2 where F 2 is a sum of degree-two functions, while 
terms of the form (3 x F\ where F\ is a sum of logarithms have to be dropped because 
they spoil the correct collinear behaviour. We also note that these correction terms are 
essential to restore the symmetry of the remainder izf under permutations of (u,v,w). 
We have checked that the two formulae for the remainder match exactly up to a factor of 
four@ 



We believe that this is the first example where the principle of maximal transcen- 



18 Note that arc the coefficients of the N 2 terms in the finite remainder defined in (5.20) of [T2"] . 
The QCD result with external gluons does not have subleading in N terms, but does have contributions 
from quarks which are proportional to powers of the number of quark flavours Np. In order to compare 
to the M = 4 form factor calculation we have suppressed such Np- dependent terms. Diagrams in the 
QCD calculation where gluons run in the loops are of course in common with the Af — 4 form factor 

(2) 

calculation, and are accounted for by A n . 

19 In this respect we note that in J\f = 4 SYM the two form factors with MHV (g + g + g~) and maximally 
non-MHV (g + g + g + ) helicity configurations are actually related, a fact that was explained in Section 3.2.3 
«f £3- 

20 The relative factors of four and two for two-loop and one-loop quantities respectively are due to 
different normalisation conventions of the gauge coupling constant. 
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dentality, first observed in [37J for anomalous dimensions of operators, relates physical 
quantities with non-trivial kinematic dependence in planar QCD and M = 4 SYM. It 
would be exciting to identify more such examples and we are confident that the technol- 
ogy of symbols will turn into a powerful tool also in QCD. 
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A Complete one-loop amplitudes with fundamental 
and adjoint generators 

At one loop, complete (planar plus non-planar) amplitudes can be written as [H] 

A^{l,...,n) = A?(l,...,n) + A$l(l,...,n) 

= N Tr(T a ^...T^)A ( ^(a 1 ,...,a n ) (A.l) 

CTGSn /2m 

|n/2j+l 

+ Yl Tr ( Ta<T1 • • • Taac ~ x )Tr(T aCTc ■ ■ • T aCT - ) AQ. (<7i, . . . , a n ) . 

0~£S n /Sn;c C=2 

Here S n is the set of permutations of n objects, and S n - C is the subset of permutations 
which leaves the double-trace structure in (lA.lj) invariant. A^\(ai, . . . ,a n ) are colour- 
ordered one-loop amplitudes, whereas A^c(ci, • • • , cr n ) are certain linear combinations 
thereof, constructed as [2] 

4Vc(l, 2, ■ ■ ■ , c - 1; c, c + 1, . . . , n) = (-1)^ 1 £ A%{a) , (A.2) 

(reCOP{a}{/3} 

where a« G {a} = {c — 1, c — 2, . . . , 2, 1}, /3j G {f3} = {c, c + 1, . . . , n — 1, n}, and 
COP{a}{/3} denotes the set of all permutations of {1,2, ... ,n} where n is held fixed, 
and such that that the cyclic ordering of the «j within {a} and of the {f3i} within is 
maintained, while allowing for all possible relative orderings of the with respect to the 

A- 
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One can also work with the alternative, remarkably compact representation of the 
complete one-loop amplitudes presented in |46j . 

A^(l,...,n) = Yl Tr(F^---F^)AW(* u ...,v n ), (A.3) 

where Fg c are adjoint generators^ The sum in (I A . 3 [) contains (n — l)!/2 terms. Note 
that planar and non-planar contributions are both contained in (1A.3[) . where only colour- 
ordered amplitudes are summed. 
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